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SUMMARY 


Experimental data are presented on the effect of Reynolds number on unsteady pressures induced by the 
pitching motion of an oscillating airfoil. Scale effects are discussed with reference to a conventional 
airfoil (NACA 64A010) and a supercritical airfoil (NLR 7301) at mean-flow conditions that support both weak 
and strong shock waves. During the experiment the Reynolds number was varied from 3 « 10 b to 12 « 10 6 at 
a Mach number and incidence necessary to Induce the required flow. Both fundamental frequency and complete 
time history data are presented over the range of reduced frequencies that Is Important in aeroelastlc 
applications. The experimental data show that viscous effects arc Important in the case of the supercrit- 
ical airfoil at all flow conditions and in the case of the conventional airfoil under strong shock-wave 
conditions. Some frequency-dependent viscous effects were also observed. 

LIST OF SYMBOLS 


C 

C 

C 


(x.t) 


P.o 


(X) 


c 

exp(-iut) 

f 


k 

M 

P(x.t) 

PINF 

QINf 


static pressure coefficient, 

(P - PINF ) /QINF 

Instantaneous pressure coefficient 

first harmonic complex amplitude of 
the unsteady pressure, per radian 

chord of wing, 0.5 m 

cos ut - 1 sin wt 

frequency, Hz, fT * 1 

qth moment of the instantaneous 
pressure coefficient 

first harmonic complex amplitude of 
the upper surface loading, per radian 

reduced frequency, uc/2U 

free-stream Mach number 

instantaneous pressure, N/m 2 

free-stream static pressure, N/m 2 

free-stream dynamic pressure, N/m 2 


Re chord Reynolds number 

T period of the motion, sec 

t time, sec 

U free-stream velocity, m/sec 

x distance along airfoil, m 

a complex amplitude of the unsteady angle 

of attack 

a mean anqle of attack 

m 3 

n( t ) instantaneous angle of attack 

u radian frequency, 1/sec 

Complex notations : 

Iml ) imaginary part of [ 1 

Magi 1 magnitude of I ) 

Ph[ ) phase of 1 1, deg 

Re( 1 real part of ( ) 


1. INTRODUCTION 

Scale effects have been considered an important element in aerodynamics research for many years. In 
early experimental studies of steady transonic flows in the 1940s, the Reynolds number, as it affected the 
state of the boundary layer approaching a shock wave, was recognized as a critical parameter (Ref. 1). In 
the 1950s the boundary-layer-trip technique was developed to simulate high Reynolds number flows in the 
wind tunnel (Ref. 2). With the development of thicker and more highly loaded airfoil sections in the 1960s, 
more complicated scale effects due to local flow separations were discovered (Ref. 3), and extensive experi- 
mental programs were developed to quantify these viscous Interactions (Refs. 4, 5). Even today, a con- 
certed effort is under way to understand the effect of Reynolds number on modern supercritical airfoil 
sections (see Gessow's introductory remarks to a recent NASA conference. Ref. 6). 

In the field of unsteady transonic aerodynamics, the scale effect has hardly been seriously considered, 
either experimentally or analytically. In the original chapter on boundary- layer effects in the AGARO 
Manual on Aeroelasticlty (Ref. 7), the authors of that work found no studies that considered the coupling 
of a boundary layer to the unsteady-pressure field. In Jones 1 review of unsteady aerodynamics in 1963 
(Ref. 8), the need for such research was reiterated; current review articles continue to cite the need for 
such research (Refs. 9. 10). 

Recently, numerical solutions to the Navier-Stokes equations were applied to unsteady transonic flow 
problems (Refs. 11-13). lhey were used to model passively excited oscillations where good qualitative 
agreement with experiment was demonstrated. Computations for a forced oscillation problem, including 
viscous effects, were reported In Ref. 14. All these codes, mostly concerned with mild transonic inter- 
actions where simple turbulence models suffice, are limited by long executing times, even on powerful com- 
puters. Future applications of the numerical method, guided by experimental data, better turbulence models. 
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and the availability of faster computers will certainly contribute to our understanding of scale effects on 
unsteady aerodynamic response. 

Experimentally, Tljdeman (Ref. 15) analyzed some of the NLR oscillatory data for scale effect. He 
accounted for the boundary layer by including the displacement effect In the calculated pressure distri- 
butions. For subsonic flows, his results show that thickness and boundary-layer effects tend to counteract 
one another In the sense that measured unsteady-pressure distributions tend to follow classical flat-plate 
theory. For transonic flows, his analysis showed that the boundary layer had a profound effect on both 
the mean shock-wave location and the characteristic unsteady-pressure peak at the mean shock locations. 
Tljdeman showed that Invlscld theories are Inadequate for predicting unsteady transonic flows, but his 
data, measured in an atmospheric wind tunnel, did not address the effect of Reynolds number directly. 

Similar comparisons and conclusions were reported by Grenon et al. (Ref. 16) In their oscillating flap 
studies. 

In this paper experimental data from tests In the 11- by 11-Foot Transonic Wind Tunnel at Ames 
Research Center are used to Illustrate the Importance of scale effects on the unsteady aerodynamics of 
both conventional and supercritical airfoils. First, the integrated upper-surface unsteady-pressure dis- 
tributions are used to assess global variations (Sec. 3.1). More detafled Information is presented with 
the aid of the fundamental frequency pressure data (Sec. 3.2), and ultimately with the Instantaneous time 
histories (Sec. 3.3). Many of the effects Involve complex unsteady viscous Interactions that are not yet 
completely understood nor easily modeled; as a result, theoretical comparisons will be restricted to 
linearized, Invlscid, unsteady aerodynamic theory. 

2. EXPERIMENTAL APPARATUS AND TEST CONDITIONS 

The experiment was conducted during the early part of 1978. The 11- by 11-Foot Transonic Wind Tunnel 
is a continuous flow facility that can be pressurized between 50 and 200 kPa for Independent control of Mach 
and Reynolds number. The test wings - an NACA 64A010 conventional airfoil section and an NLR 7301 super- 
critical airfoil section - were mounted between two floor-to-celling splitter plates Installed in the 
3.35- by 3.35-m slotted test section. The 0.5-m-chord by 1.35-m-span wings were fabricated from a light- 
weight graphite-epoxy composite material. A sketch of the test apparatus is shown in Fig. 1, and a complete 
description of the test hardware, motion generators, and models is given in Ref. 17. The chord Reynolds 
number range was approximately 3 * 10 6 to 12 * 10 s . No boundary- layer trips were used in this test. 

The data reported in this paper form a small portion of the extensive data base that was collected 
during the test. Some of the data were previously reported in Refs. 18 and 19, and a discussion of the 
method used to acquire and validate the unsteady data is presented in Ref. 20. The data subset relating to 
scale effects is presented in Table 1. Two broad categories are represented: conditions 1 to 3 represent 

attached flows and conditions 4 and 5 are examples of separated flows. As will be shown presently, both 
classes of flows can exhibit scale effects, with the separated flow data indicating a more erratic behavior. 

3. PRESENTATION OF DATA AND DISCUSSION 

3.1 Variation of the Global Parameter I 0 , a With Frequency 

A convenient measure of the load on the upper surface is the chordwise integral of the unsteady- 
pressure distribution. As introduced in Ref. 19, this parameter is defined by 

I (t) ■= / C (x,t)(x/c) q d(x/c) 
q o p 

where q = 0 is a measure of the lift attributed to the upper surface (with sign reversed). A Fourier 
decomposition of I 0 (t) defines the amplitude of the first harmonic component 

aI o.a = f / I o ( t ) exp ( - i u>t ) dt 
0 

Considering first the attached flows, conditions 1-3 in Table 1, the variation of the magnitude of I 0 , a 
with frequency is showr in Table 2, with Reynolds number as parameter. Values from subsonic, inviscid, 
unsteady aerodynamic theory (Ref. 21) are also shown for comparison. The frequency variation, typical of 
attached transonic flows - decreasing amplitude with increasing frequency - persists at all Reynolds num- 
bers. The data also show decreasing amplitudes with increasing Reynolds numbers at most frequencies. For 
the subsonic flow (M = 0.500) the variations are minimal and can probably be accounted for by simple inclu- 
sion of the boundary-layer displacement. The transonic flows, being sensitive to shock-wave/boundary-layer 
interactions, show slightly more severe scale effects. 

A graphic summary of the scale effect for conditions 1-3 is shown in Fig. 2 where the variation of the 
complex amplitude with Reynolds number is presented at a reduced frequency k = 0.05. As expected, the 
subsonic flow condition is least affected by Reynolds number. The trend for both the conventional and 
supercritical data at transonic Mach numbers is toward decreasing in-phase and out-of-phase (Re and Im) 
components with Reynolds number. It is noteworthy that, for each flow condition, the phase does not change 
by more than 5’ over the Reynolds numbers range considered. 

In Refs. 18 and 19 it was shown that the aerodynamic transfer function (variation of load with fre- 
quency) was not even qualitatively consistent with linear theory when the shock wave was strong enough to 
separate the boundary layer. The separated boundary layer had a major effect on the loads, and it is 
expected that scale effects would also be very important. Table 3 shows the magnitude of the complex 
amplitudes for conditions where strong unsteady shock-wave/boundary-layer interactions are encountered. 

No discernible trends are evident from these data, in fact, the magnitudes undergo alarming variations 




with Reynolds number when compared with the attached-flow data. Due to this complex behavior, further 
information must be gleaned from the unsteady-pressure data Itself. 

3.2 Unsteady-Pressure Distributions 
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The first harmonic complex unsteady pressure Is defined by the following expression: 

aC (x) * i f C (x,t)exp(-1wt)dt 
P.“ T » 0 P 

An analysis of the In-phase and out-of-phase (Re and Im) component of Cp 1? will show that the distribu- 
tions for attached-flow conditions are qualitatively consistent with the simplest theoretical models that 
include viscous effects by consideration of the boundary- layer displacement. The separated flow, however, 
shows important unsteady effects that are attributed to both the steady-flow scale effect and the unsteady 
viscous response of the boundary layer. 

The data for attached-flow conditions 1-3 are shown in Figs. 3-5, respectively. For each frequency, 
the mean, in-phase, and out-of-phase components of the unsteady-pressure distribution are presented. 

Figure 3 shows that the subsonic flow (o = 0, M = 0.5) has minimal scale effect. The only difference 
between Reynolds numbers is the slight dip in the in-phase component at wing station x/c * 0.5. Subli- 
mation photographs showed that transition occurred at this location at Re « 2.5 * 10 s ; at Re = 10 * 10 s , 
transition was very close to the leading edge. In any event, the state of the boundary layer has no effect 
on the unsteady-pressure distribution, which assumes the classical shape that was predicted many years ago 
from linearized subsonic theory. The same airfoil at transonic speeds in Fig. 4 (a = 0, M = 0.8) shows 
that scale effect is also quite minimal, Aside from a slight upwind movement of the mean shock position 
that affects the chordwise location of the ui. steady-pressure peak, the effect is minor. Sublimation photo- 
graphs at this condition showed that transition occurs at the shock wave ut Re = 3.3 x 10 6 ; while leading- 
edge transition was observed at Re * 12.6 * 10 6 . The last attached flow condition, shown in Fig. 5 for 
the NLR 7301 supercritical airfoil, exhibits mo-e severe scale effects. Although t' e region of rapid com- 
pression (probably not a shock wave at this supercritical design condition) moves only slightly, the 
unsteady pressure is quite different at the three Reynolds numbers indicated. These differences me prob- 
ably due to the large extent of supercritical flow on the upper surface that is affected by tne change in 
effective airfoil shape due to the unsteady boundary- layer growth. The difficulty in distinguishing dif- 
ferences by examining the integrated values shown previously is obvious because of the varying contribution 
to the loads from positive and negative lobes of the unsteady pressures. 

The data for shock-induced separation are more difficult to analyze. These data are shown in Figs. 6 
and 7. Figure 6 shows the large effect of Reynolds number for the NACA 64A010 at a - 4®, M = 0.8 (condi- 
tion 4 in Table 1). The mean flow differs not so much in the shock position as in the extent of separated 
flow downstream of the shock. This separated flow has a marked effect on the in-phase unsteady pressure at 
the shock wave and beyond. The out-of-phase pressures are not very much different at the two Reynolds num- 
bers srown. The mechanisms whereby the in-phase pressures are more sensitive than the out-of-phase pres- 
sures have yet to be explained. Figure 7 shows similar data for the supercritical airfoil at a Mach number 

beyond its design point (condition 5 in Table 1). The mean flow differs from the previous case in the 

absence of any discernible training-edge pressure recovery downstream of the shock wave. This indicates 
a more severe flow separation. The data for the unsteady components of the pressure distribution also show 
a significant sensitivity to Reynolds number. For the in-phase component, especially, variations in the 
negative contribution to the load change rapidly with Reynolds number. This large variation gives rise to 
the confusing trends in the unsteady loads shown in Table 3. Again, it should be noted that the in-phase 
components appear to be more sensitive to Reynolds number than the out-of-phase components. 

It is clear from the preceding discussion that a more fundamental data set needs to be examined to 
clarify the underlying physical mechanisms behind the scale effect. The final sequence of data will show 

how the Reynolds number effects the instantaneous unsteady-pressure time histories. 

3.3 Instantaneous Pressure Data 

The measured unsteady-pressure coefficients Cn(x,t) are presented in Figs. 8 to 12 for the five flow 
conditions listed in Table 1. At each chordwise station (identified by a numerical key), data are shown at 
two Reynolds numbers. The mean portion of the instantaneous pressure coefficient was suppressed for clarity. 
(The reference line for each trace is the corresponding tic mark on the airfoil contour.) As explained in 
Ref. 23, the dynamic data have been processed to eliminate all asynchronous signals. Thus, all of the dips 
and bulges shown in the data are truly periodic. Whether they can be traced to a particular fluid- 

mechanical event at that Instant is argumentative, but they are included here for completeness. The input 

motion is the same for both Reynolds numbers and can be used as the phase reference (wt = 0 when o • » maJ( ). 

For the most part, data are shown at the low, reduced frequency k * 0.05. 

The subsonic flow (condition 1) is presented in Fig. 8, and, as expected, scale effect is minimal. 

The usual trend of decreasing amplitude with increasing chordwise location is obvious. The distorted sig- 
nal at Re * 2.5 « 10 6 in traces 9 and 10 is apparently caused by the transitional boundary layer. The 
"noise" on traces 15 and 16 at low Reynolds number has no discernible fluid-mechanical origin. 

The effect of increasing the Mach number to 0.80 is shown in Fig. 9. The presence of the shock wave 

is apparent by its distortion of the pressure signal of both Reynolds numbers, causing a severe local scale 

effect. However, both upstream and downstream of the shock wave, the Reynolds number does not have a sig- 
nificant effect. The global ramifications of these local effects were tabulated in Table 2. 

Data from the supercritical airfoil at its shock-free design condition are shown in Figs. 10a and 10b. 

Data are presented at two frequencies to show the complicated cross-coupling between frequency and scale 

effects that was not present in the previous cases. In Fig. 10a data at the low reduced frequency of 
k =■ 0 05 show very significant scale effejts in the supercritical flow region. At R<; * 12.6 x 10 s there 
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Is a higher amplitude, more distortion, and significant phase shifting. Whether these effects ctn be 
explained by simple boundary- layer displacement corrections renal ns to be answered. In Fig. 10b there Is 
« surprising trend toward reduced harmonic distortions at higher Reynolds numbers (traces 5-9), although 
there are still significant phase shifts attributable to scale effect. At the "shock wave* (which appears 
during the cycle), the effect of Reynolds number 1$ apparent. Consistent with the previous attached flow 
data, there Is minimal response near the trailing edge at both frequencies. 

Data pertaining to the second flow regime (e.g. , shock-induced separation), are presented In Figs. 11 
and 12. In Fig. 11 the scale effect Is most prominent at the shock-wave location (traces £ and 9). The 
fundamental frequency data shown In Fig. 6 Indicate that the In-phase component Is most severely affected 
by the Reynolds number. The large change In the unsteady pressure Is apparent. Upstream, there Is a sig- 
nificant Increased phase lag with Increasing Reynolds number. This phase lag persists Into the separated- 
flow region aft of the shock. Unfortunately, no data are available at k ■ 0.05. 

In Fig. 12 the well -separated flow over the supercritical airfoil does not Induce large unsteady dis- 
turbances downstream of the shock. This seems to ameliorate the scale effect in the leading-edge region, 
but not at the shock itself. Data at k * 0.05 show similar trends. 

This cursory examination of the unsteady-pressure traces Indicates the sensitivity of the flow patterns 
to geometry, frequency, mean-flow conditions, and Reynolds numoer. It Is clear that predictive schemes must 
Include, as a minimum, physical models of all of these parameters. 

4. CONCLUDING REMARKS 

Examination of some of the data from a series of tests on oscillating airfoils In the Ames 11- by 11- 
Foot Transonic Wind Tunnel showed that scale effects could be a very Important factor In the unsteady aero- 
dynamic behavior of the airfoils. Configurations with mean flows that support fully attached boundary 
layers and weak shock waves have unsteady response characteristics that are mild functions of Reynolds 
number. These configurations can probably be modeled with currently available computational tools. An 
exception, perhaps. Is airfoils that possess mean flows with extensive regions of supercritical flow. The 
detailed pressure distributions contain substantial scale effects, but the overall loads may not be so 
severely affected. 

Configurations with mean flows having detached boundary layers are enigmatic. As reported In Refs. 18 
and 19, the contribution to the unsteady load may be caused by a delicate balancing of positive and nega- 
tive lobes in the unsteady pressures (see Fig. 7). For these conditions, shape changes caused by scale 
effects may have severe ramifications. The erratic behavior of the unsteady loads shown In Table 3 Is 
Illustrative of this effect. The major technological application of unsteady aerodynamics Is to aeroelastic 
analysis and design. The sensitivity of flutter boundaries and stability margins to scale effects must be 
considered in those situations. 

In the mathematical modeling of unsteady transonic aerodynamics, It Is clear that the Reynolds number 
needs to be included as a primary parameter. Progress In computational research during the past decade was 
characterized by the inclusion of nonlinear thickness effects, and It Is hoped that this and other experi- 
mental data, along with new computational efforts, will advance our knowledge well beyond the Invlsdd 
approximation. 
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TABLE 1.- RANGE OF FLOW CONDITIONS CONSIDERED 


Flow condition 

Airfoil 

M 


Motion 

1 

NACA 64A010 

0.500 

0 

Pitching tl* at 0,25 c 

2 

NACA 64A010 

0.796 

0 

Pitching *1* at 0.25 c 

3 

NLR 7301 

0.752 

0.37 

Pitching ±0.5' at C.a.j c 

4 

NACA 64 AO 10 

0.789 

4.0 

Pitching at 0.25 c 

5 

NLR 7301 

0.807 

0.38 

Pitting cQ.5“ at C.40 c 


TABLE 2.- MAGNITUDE OF UNSTEADY LOADING ATTRIBUTED TO UPPER SURFACE - ATTACHED FLOW 


Re « 10-6 




Mag 





k * 0 

k * 0.025 k - 0.05 k 

• 0.10 

k « 0.15 

k « 0.20 

k * 0.25 

k * 0.30 




Flow condition 1 (M * 

■ 0.500) 




2.5 

3.7 

3.67 

3.59 

3.27 

2.97 

2.83 


_ 

5 

3.8 

- 

- 

- 

- 

2.45 

- 

- 

10 

- 

- 

3.43 

3.19 

2.91 

2.74 

2.58 

- 

Inviscid theory 

3.64 

3.42 

3.22 

2.91 

2.69 

2.52 

2.44 

2.38 




Flow conditior 

i 2 (M - 

0.796) 




3.3 

5.5 

_ 

4.87 

_ 

3.19 

2.84 

2.24 


6.7 

5.7 

- 

- 

- 

- 

2.67 

- 

- 

12.6 

- 

4.62 

4.41 

3.80 

2.90 

2.68 

2.16 

2.16 

Inviscid theory 

5.25 

4.61 

4.12 

3.46 

3.07 

2.84 

2.71 

2.63 




Flow condition 3 (M * 

0.752) 




3.3 

- 

- 

6.46 

. 

_ 

2.98 



6.2 

- 

6.56 

6.19 

4.80 

- 

3.67 

- 

2.01 

11.5 

- 

6.05 

5.62 

4.63 

3.65 

3.12 

- 

1.94 

Inviscid theory 

4.76 

4.28 

3.88 

3.22 

'.97 

2.75 

2.62 

2.54 


TABLE 3.- MAGNITUDE OF UNSTEADY LOADING ATTRIBUTED TO UPPER SURF.\: C - SEPARATED FLOW 


Mag Ij.j 

Re , 10-6 — — - 

k = 0.05 k = 0.20 


Flow condition 4 (M * 0.789) 


6.2 



4.48 

11.9 


1.67 

5.57 

Inviscid 

theory 

4.12 

2.64 


Flow condition 

5 (M -- 

0.80/,, 

3.3 


2.39 

1.53 

6.3 


0.44 

0.69 

11.7 


1.08 

0.88 

Inviscid 

theory 

4.13 

2.82 
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40 STATIC PRESSURE 



IINPUTI 


Fig. 1. Oscillating airfoil test apparatus Installed 
In the 11- by 11-Foot (3.4 by 3.4 m) Transonic 
Mind Tunnel at Ames Research Center. 
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Fig. 2. Variation of the unsteady loading attributed 
to the upper surface with Reynolds number; 
low-frequency data, k * 0.05. 
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